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Starting with Wannier's' and Anderson's^ seminal papers 
the triangular-lattice antiferromagnet has attracted a lot of 
attention. One spectacular feature of this model system for 
a strongly frustrated magnet is the unconventional magne- 
tization process of the triangular-lattice Heisenberg antifer- 
romagnet.''"'^ While for the classical isotropic Heisenberg 
model at zero temperature the magnetization M increases lin- 
early with the applied magnetic field H, thermal or quantum 
fluctuations induce a plateau at 1/3 of the saturation mag- 
netization Msat ('order from disorder' phenomenon). For the 
extreme quantum case, i.e. spin quantum number 5=1/2, 
this plateau at T = has been widely discussed. Very re- 
cently an almost perfect experimental realization of an 5 = 1 /2 
triangular-lattice Heisenberg antiferromagnet has been re- 
ported for Ba3CoSb209,'^ where the entire measured magne- 
tization curve M(H) including the 1/3 plateau is in excellent 
agreement with theoretical predictions. In particular, the the- 
oretical magnetization data obtained by high-order coupled- 
cluster approximation" and by large-scale exact diagonaliza- 
tion approach'^ almost coincide with the measured ones over 
a wide range of the magnetic field. 

Another interesting triangular-lattice magnet is 
Ba3NiSb209 that is considered as a good candidate of 
an i= 1 triangular-lattice Heisenberg antiferromagnet.'^' 
In Ba3NiSb209, the uniform triangular lattice of magnetic 
ions is realized as in Ba3CoSb209. Indeed, recently it 
has been found that the magnetization curve of this com- 
pound also exhibits a clear 1/3 plateau.'^ By contrast to 
the well-investigated spin-half case there are much less 
reliable theoretical studies of the s = \ model relevant for 
Ba3NiSb209. Motivated by the excellent agreement of 
theoretical predictions based on exact diagonalization (ED) 
and coupled-cluster method (CCM) and experimental results 
for the spin-half compound Ba3CoSb209 reported in Ref. 17 
we apply in this paper these methods to the s- \ model to 
provide theoretical data to compare with experimental results. 

For the s- \ model the Lanczos ED for = 27 sites and the 
CCM-SUB«-« approximation-" for « = 2,4, 6,8 are used to 
calculate the field-dependent properties of the model. In this 
short note we will not explain details of the Lanczos ED and 
the CCM. We refer the interested reader e.g. to Refs. 15,21,22 
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Fig. 1. The magnetization curve of the s= \ triangular-lattice Heisenberg 
antiferromagnet calculated by the CCM and by ED compared with the ex- 
perimental data for Ba3NiSB209 measured at T = 1.3 K, see Ref. 16. For the 
compaiison we have plotted the theoretical data using the exchange parame- 
ter Jjk^ = 21.9 K and the g factor g = 2.392. The inset shows magnetization 
curves up to the saturation calculated by the CCM and by ED. 
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Fig. 2. The derivative susceptibility dM/dH of the s=i triangular-lattice 
Heisenberg antifeiTomagnet calculated by the CCM compared with the ex- 
perimental data for Ba3NiSB209 measured at T = 1.3 K, see Ref. 16. For the 
comparison we have plotted the theoretical data using the exchange parame- 
ter J/ke = 21.9 K and the g factor g = 2.392. 



and 1 1, 23-25, respectively. 

First we report the main theoretical results. The satura- 
tion field is //sat = 9Js = 97, where J is the nearest-neighbor 
exchange coupling. The zero-field uniform susceptibility x 
calculated with CCM-SUB«-« approximation is = 0.10790, 
0.09932, 0.09785, and 0.09679 for «=2, 4, 6, and 8, re- 
spectively. Moreover, it is useful to extrapolate the "raw" 
SUB«-« data to «— >oo by ;if(«)=co + ci(l/n) + C2(l/n)^ 
which yields finally our CCM estimate for the suscep- 
tibility xccM - 0.0956. This value is in exceUent agree- 
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Fig. 3. Upper panel: Spin orderings in magnetic field. Lower panel: CCM- 
SUB6-6 data for the sublattice magnetizations Ma,b(H) and Mc{H) (left) and 
the canting angles a(H)/7T and /3(H)/7r (right) as functions of the magnetic 
field H. 



ment with the spin-wave resuh, ° ;^'sw =0.0953. (We re- 
mark that this value of xs in Ref. 26 was referred to 
as x± in this article and furthermore that it was defined 
per volume. The left and right endpoints of the 1/3- 
plateau, H^i and H^i are determined to //d/y = 2.817 and 
//c2// = 3.695 (ED, N^27), //ei/^ = 2.727, 2.788, 2.809, 
2.8142 and H^2/J = 3.617, 3.674, 3.648, 3.637 (CCM, SUB«- 
n for n-2, 4, 6, 8, respectively). As for the susceptibil- 
ity we may extrapolate to n^oo yielding //ci/-/ = 2.839 
and Hc2/J-3.552. We can compare these values with 
corresponding ones obtained by linear spin-wave theory,"^ 
//ci/y = 2.748 and Hc2/J - 3.645. Improving the linear spin- 
wave approximation by using a self-consistent approach 
Takano et al.^^ obtained //ci// = 2.760 and //^a/^ = 3.585, 
which are closer to our CCM estimates. The sublattice mag- 
netization in the collinear plateau 'up-up-down' state (spin 
configuration II in Fig. 3) calculated within CCM-SUB«-« 
approximation are M^p =Ma.b =0.93251, 0.91565, 0.91188, 
0.91097 and Mjown = Mc = - 0.86502, - 0.83129, - 0.82375, 
-0.82193 for n = 2,4, 6,8, respectively. Again we may ex- 
trapolate to « — > oo yielding M^p - 0.909 and Mdown = - 0.817. 

Similar to Ref. 17 we now compare directly the theoret- 
ical predictions for the ground state with low-temperature 
experimental results reported in Ref. 16, see Figs. 1 and 2, 
where experimental data were corrected for the Van Vleck 
paramagnetic susceptibility of xvv ~ 2.4 x 10"^ emu/mol = 
4.3 X 10 '•//B/(Ni^+ T), which has not been done in Ref. 16. 
Due to the finite-temperature effect, small single-ion 
anisotropy and small anisotropy of the g factor the experi- 
mental curve is smeared out around the endpoints of the 1/3- 
plateau, Hd and H^i. Nonetheless, the excellent agreement 
between theory and experiment is obvious. A good agreement 
between experimental and theoretical data is also obtained 
for the derivative susceptibility dM/dH, shown in Fig. 2. 
We may conclude, that the presented data demonstrate that 
Ba3NiSb2099 is well described by the s-l triangular-lattice 



Heisenberg antiferromagnet. 

In addition to the total magnetization our CCM approach 

allows to calculate the sublattice magnetizations as well as the 
canting angles (see Fig. 3). The three different ground-state 
regimes illustrated in the upper panel of Fig. 3 are clearly seen 
in the behavior of these quantities. While the total magneti- 
zation M = 2Mup + Mdown increases monotonously with mag- 
netic field (except within the plateau region), both sublattice 
magnetizations, Myp and Mdown, exhibit regions of negative 
slope dMsh/dH <0. 
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